Heavy quark diffusion in perturbative QCD at next-to-leading order 
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We compute the momentum diffusion coefficient of a nonrelativistic heavy quark in a hot QCD 
plasma, to next-to-leading order in the weak coupling expansion. Corrections arise at 0{gs); phys- 
ically they represent interference between overlapping scatterings, as well as soft, electric scale 
(p ^ gT) gauge field physics, which we treat using the hard thermal loop (HTL) effective theory. 
In 3-color, 3-flavor QCD, the momentum diffusion constant of a fundamental representation heavy 
quark at NLO is k = i|^a?T^(ln^ + 0.07428 + 1.8869ffs). The convergence of the weak coupling 
expansion is poor. 
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The experimental program at RHIC and the future 
heavy ion program at the LHC are exploring the behav- 
ior of the QCD plasma at temperatures above the "de- 
confinement" temperature of ^ 170 MeV. So far the evi- 
dence is for a medium which interacts more strongly and 
thermalizes more quickly than expected. For instance, 
experimental results on elliptic flow are well explained 
by hydrodynamics [1] but only if the shear viscosity is 
much less than a naive extrapolation of weak coupling 
calculations [2, 3]. Similarly, heavy quarks display sub- 
stantial elliptic flow and a degraded energy spectrum [4] , 
implying stronger medium interactions than extrapolated 
weak-coupling calculations can easily accommodate [5]. 

This raises the general question; how well can we trust 
weak coupling calculations for dynamical quantities in 
hot QCD at couplings anywhere close to those probed in 
experiments? Naively the perturbative series converges 
in (possibly non-integer) powers of the strong coupling, 
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Q!s = 1^ ~ 0.4 for relevant temperatures. But perturba- 
tive series often show convergence which is much better 
or much worse than one would guess from the value of 
the coupling. In general determining how well a per- 
turbative expansion converges requires evaluating a few 
terms in the expansion. Unfortunately, no dynamical 
transport quantity in QCD which involves large length 
or time scales (such as shear and bulk viscosity, elec- 
tric conductivity, photon production, and heavy quark 
momentum diffusion and energy loss) is known beyond 
leading order. Most of the leading-order calculations are 
recent and quite involved. 

Here we present a next-to-leading order calculation of 
the theoretically simplest of these quantities, the mo- 
mentum diffusion coefficient of a nonrelativistic heavy 
quark. This coefficient (partially) characterizes how 
quickly heavy quarks are thermalized and swept up in 
the collective flow of the plasma. 

A heavy quark, M ^ T, in or near equilibrium has a 
typical momentum p ~ V AIT ^ T large compared to 
the plasma scale and it therefore takes a parametrically 
long time for the momentum to change appreciably. This 
means that momentum changes accumulate from many 



uncorrelated "kicks," so on long time scales p will evolve 
via Langevin dynamics. 
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The relaxation rate rjjj and the momentum diffusion con- 
stant K are related by a fluctuation-dissipation relation, 
rjo — 2MT ' '^hich follows on general thermodynami- 
cal grounds. Thus the dynamics of the nonrelativistic 
heavy quark is completely set by the single parameter n. 
This parameter can be obtained by computing the mean 
squared momentum transfer per unit time in the under- 
lying microscopic theory. In gauge theory, this mean 
squared momentum transfer equals the time integrated 
correlator of two electric field operators connected by fun- 
damental Wilson lines [6] : 

n={- dTTTH{W{T;0)^ EfiT)t%WiT;0)E'liO)t'ij), 
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where W{T; 0) denotes a fundamental Wilson line run- 
ning from t — to T along the static trajectory of the 
heavy quark and = 3 is the dimension of the heavy 
quark's representation. 

Intuitively, Eq. (2) is exactly the force-force correlator 
of Eq. (1), with the forces given by electric fields and 
the Wilson line representing the gauge rotation of the 
heavy quark due to propagation, which ensures gauge 
invariance. Because of operator ordering issues, the Wil- 
son lines shown are not equivalent to connecting the E 
fields with an adjoint Wilson line, and in fact such Wil- 
son lines are even required in QED (diffusion of ions in 
a QED plasma depends on the ionic charge Z in a more 
complicated way than Z'^ only because of these Wilson 
lines, which account for the reaction of the plasma to the 
presence of the charge) . 

We start by showing how this formula reproduces the 
well known [7] leading order momentum diffusion coef- 
ficient. At this order, (2) simplifies to a zero-frequency 
Wightman correlator of two A° fields (the fields do not 
contribute to the electric field operators at zero- frequency 
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FIG. 1: Leading-order contribution to heavy quark diffusion 

and its correspondence to scattering processes. On the left 
the double line represents the Wilson line; on the right it is 
the heavy quark external states. 



in covariant and Coulomb gauges [14]) 
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where Ch = | is the Casimir of the heavy quark's repre- 
sentation. This Wightman correlator can be evaluated in 
terms of the squared matrix elements of t-channel scat- 
tering processes involving the heavy quark, as illustrated 
in Fig. 1. These are the only processes which contribute 
in our case, Compton-like processes being suppressed in 
the low velocity limit. The result reduces to [8] 



K = 
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dp 
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Here p is the transferred momentum and q is the energy 
of the light scattering target. Since the heavy quark is at 
rest, the initial and final light-particle energies arc equal 
and p is purely spatial, which is why the medium modifi- 
cation of the exchanged gluon propagator is purely Dcbyc 
screening with a Debye mass = g'^T'^{Nc + Nf/2)/3. 
The inclusion of these HTL corrections is essential for 
obtaining the complete leading order result, otherwise 
K would be infrared divergent in the region of soft mo- 
mentum transfer p. Formally taking niu ^ T, the in- 
tegral is dominated by g ~ T and p in the parametric 
range nio <^ p <, T. The strict leading-order evaluation 
of Eq. (4) yields 
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-0.64718. 



with e = i - 7. + ^ 

When the exchange momentum p is hard, p > T, then 
higher loop corrections to the propagators and vertices in 
Fig. 1 represent 0{g^) corrections. However, the expres- 
sion (4) for At receives an 0{g) contribution from scat- 
terings against soft gluons, q ~ mo ■ Both the dispersion 
relations and the interactions of such gluons are modified 
at the 0(1) level; at leading order these modifications are 



described by hard thermal loops. Therefore there will be 
0{g) corrections to the above calculation. But this is 
not the only source of 0{g) next-to-leading order (NLO) 
corrections. 

Another source is associated with overlapping scatter- 
ing events: the total scattering rate for a hard particle is 
~ 5^r, and is dominated by t-Channel Coulombic scat- 
terings involving soft momentum transfers. These soft 
scatterings have a duration of order 1 / mo 1 / gT and 
therefore there is an 0{g) probability that two such scat- 
tering events overlap with each other. This is relevant 
in QCD (though not in QED, see below) because each 
scattering color-rotates the participants. 
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FIG. 2: Diagrams required at NLO. The double line is the 
Wilson line; otherwise all propagators are soft and HTL re- 
summed and all vertices include the HTL vertex. All lines 
attached to the Wilson line are longitudinal. 

We need a systematic way of evaluating these NLO 
effects. This is provided by a loopwise expansion for 
Eq. (2). The diagrams needed at NLO are shown in 
Fig. 2. The diagrammatic series is convergent in pow- 
ers of g provided one incorporates HTL corrections in 
propagators and vertices wherever momenta are soft [9], 
unless a diagram is sensitive to the magnetic scale ^ g^T, 
which would be signaled by an infrared divergence in the 
evaluation of a Fcynman diagram. This does not occur in 
the current calculation; the diagrams shown in Fig. 2 are 
all IR and UV convergent, after the leading-order con- 
tribution is subtracted off from the transverse, pole-pole 
contribution of diagram (A). Since the momenta are soft, 
the ordering issues for the Wilson lines are subdominant 
and we may replace the two Wilson lines in Eq. (2) with 
an adjoint Wilson line; all diagrams involve the group 
theoretic combination ChCa and we may represent the 
NLO correction as the coefficient C defined by 
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with 0{g^) corrections. There is no 0{g) NLO correction 

in QED, where the (bare and HTL) vertices involved in 
diagrams (A), {B), (C) do not exist and the Wilson fine 
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in {D) is trivial. 

Since Eq. (2) involves unequal time correlators we have 
found it most convenient to evaluate it in the real-time 
(Schwingcr-Kcldysh) formalism. This required an exten- 
sion of the HTL formalism to the closed time path in the 
Schwinger-Keldysh (r, a) basis [10], which is convenient 
for treating soft physics because Bosc-Einstein factors 
only arise in one propagator. We work in strict Coulomb 
gauge. The measurable, Eq. (2), is gauge invariant and 
the HTL expansion should respect gauge invariance, so 
we expect the sum of diagrams systematically evaluated 
in powers of g to produce gauge invariant results, though 
the results for individual diagrams probably arc not. 

The effect of diagram (A) can be divided into the real 
and the imaginary part of the self-energy correction. The 
real part is the simplest to compute; it represents a cor- 
rection to the Debye mass which can actually be evalu- 
ated within the 3-D dimensionally reduced theory [11]. 
4D Coulomb gauge corresponds to 3D Landau gauge; in 
this gauge the self-energy receives a nonzero, momentum- 
dependent contribution when one propagator in the self- 
energy is transverse and the other is longitudinal {Aq in 
the 3D theory). The correction is 
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6ml = -^C^g'^T 
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The contribution to C is found by expanding {p^ + ml + 
5ml)~'^ — {p"^ + ml)'"^) ~ —25ml{p'^ + ml)~^ and finding 
the shift to Eq. (4). Straightforward integration gives 



Cr, 



(A) 
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0.77199. 



The next simplest contributions are from diagrams 
(C) and {D). Physically, (C) accounts for real and 



2 in Eq. (4), (7) C, 



virtual corrections in which the light scattcrcr under- 
goes an additional soft scattering or soft plasmon emis- 
sion/absorption. Diagram (D) is the same but for the 
heavy quark. In QED there is a cancellation between 
vertex orderings but in QCD one instead picks up a com- 
mutator of color operators. The contributions of these 
two diagrams are 

= ^^i(2.)3 il+p2)2 J^J (^^- (Q)-^--^ 

6{v ■ jP-Q)) + 6{v ■ (P+Q)) - 26{v ■ P) 

2 poo pOC f-OC 

C(D) = ^ J P'^dp J q^dq j duj 

^GO0(a;,p)-G°0(0,p) 

X Cr^^ (.w, g j . [\)) 

In writing these expressions we have scaled all mo- 
menta by m,o and scaled out all powers of T. Here 
v^^ = (l,v) and G^^ is the ordering- averaged gauge 
field correlator, related to the retarded correlator via[15] 

= (2nB(w)-fl)ReGr(^,rt - 221 RcGr(^,p). 
These expressions can be simplified somewhat but must 
be evaluated by numerical quadratures. We find [12] 



(C) 



-0.132916(1) and G(£,) = 0.067526(1). 



The most involved calculation is for the imaginary con- 
tribution of the self-energy loop in diagram {A). This 
bears some similarity to the calculation of the gluon 
damping rate by Braaten and Pisarski [13], but the "ex- 
ternal" momentum P is now spacelike. Therefore the in- 
tegrals encountered are four rather than two dimensional 
(one must integrate over p and 9pq), and the kinematics 
allow for processes involving two soft plasmons on their 
mass shells, as well as virtual corrections to the tree pro- 
cess of Fig. 1. The contribution to G can be written 
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, (10) 



where we have introduced 
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47r {vQ-ie){v R-ie) ' ^^^^ 
V = 2gV + (i? + P)^^0-(Q + P),<50 (12) 

to denote objects that enter the HTL and tree vertices. 
The evaluation is lengthy [12], and rather remarkably, 
turns out to be separable into pole-pole, pole-cut and 
cut-cut contributions, in analogy to what was found by 
Braaten and Pisarski. Two subtractions are required. 



First, as mentioned above, at large momenta the pole- 
pole contribution when both gauge bosons are transverse 
duplicates the tree process of Fig. 1; this must be sub- 
tracted. Further, evaluating the integrals in Eq. (3) for 
finite mo already incorporates an NLO correction besides 
what is in Eq. (5) . We will take the leading contribution 
to be the result including this NLO correction (which 
corresponds to Ggq. (3) = |^.) After these subtractions 
we obtain (numerically) Gii„ (^) = 0.9097(1). 

Diagram [B) involves the correlator of three Aq fields 
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connected by an HTL 3-point function (the tree vertex 
vanishes), and accounts for interference between scatter- 
ing events occurring on the Ught scatterer's side and on 
the heavy quark's side. One of the A'^ fields carries zero 
frequency, and the contributions can be organized ac- 
cording to whether the zero frequency propagator is cut 
or retarded, 

r - f d'p P f d\ M^^(q,-r) 

(B) j (2^^3 (^^^2)2 j (2^)3 (l+^2)(l+^2) Ai^j 



im (B) 



= 127r 



(Pp p^ 



(27r)3l+p2_/ (27r)4 
G?0(Q)MOO(Q,p) - GOO(Q)MOO(Q,p)* 



.(14) 



The contributions C^e (b) and Cim (b) are cfosely anal- 
ogous to the real and imaginary parts of diagram [A), 
respectively. We find [12] Cre(B) = -0.04829(1) and 
Cim(B) = -0.07338(1). 

Discussion 

The heavy quark diffusion coefficient can be computed 
beyond leading order in the weak coupling expansion. 
The first corrections arise at 0{gs) and describe "soft" 
p,uj ^ QsT physics including interference between scat- 
terings and plasma corrections to interaction strengths. 
The calculation requires the HTL effective theory. 
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FIG. 3: Comparison of leading and NLO results for iVf = 3 
QCD as a function of coupling. 

The ratio of the NLO correction to the LO result is 
independent of the representation of the heavy quark and 
is proportional to the group's adjoint Casimir Ca [Ca = 
Nc in SU(A'c) gauge theory, in QED). Numerically, we 
find (see Eq. (6)) C = 1.4946 + CEq. (3), or, for 3 flavors, 

K = -ifa^T^ (ln—+ 0.07428 -|- 1.9026 + 0{gl) I . 

(15) 

The correction is positive, meaning faster equilibration 
of heavy quarks. As shown in Fig. 3, for realistic values 



of the strong coupling the correction is large-a factor of 
2 already at as = 0.03. 

Our result suggests that the convergence of the pertur- 
bativc expansion for dynamical quantities is poor. Why? 
About ^ of the NLO coefficient in Eq. (15) (the part 
we called CEq. (3) ) is incorporated by integrating Eq. (4) 
numerically rather than expanding it into Eq. (5). An- 
other third, Cre (A) , can be approximately included into 
Eq. (4) by giving the real part of the self-energy its full p 
dependence rather than approximating it with its small 
p limit, m^. The remaining third represents complicated 
and nontrivial many-body physics. 

It would be interesting to make similar calculations for 
other transport coefficients such as shear viscosity, and 
to extend the present calculation to J\f=A SYM theory. 
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